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Abstract 

We compute the spectral index for scalar perturbations generated in a primordial 
inflationary modeh In this model, the transition of the inflationary phase to the 
radiative era is achieved through the decay of the cosmological term leading a second 
order phase transition and the characteristics of the model allow to implement a set 
of initial conditions where the perturbations display a thermal spectrum when they 
emerge from the horizon. The obtained value for the spectral index is equal to 2, a 
result that depends very weakly on the various parameters of the model and on the 
initial conditions used. 

PACS number(s): 04.20.Cv., 04.20.Me 

1 Introduction 

The big-bang standard model, based on the Friedmann-Lemaitre- Robertson- Walker (FLRW) 
solutions of Einstein's equations suffers from some drawbacks which leaded cosmolo- 
gists to propose the so-called inflationary models |0. The common element to all these 
models is an exponential expansion phase of the size of the Universe prior to the beginning 
of the standard scenario that describes the adiabatic era extending from temperatures of 
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order of 10 ^ at Grand Unification era (GU) until 10 today's temperature, in Planck's 
units (which are the units that will be used throughout all this work). 

From a kinetic point of view, the exponential expansion phase corresponds to a de 
Sitter geometry which later on will transform itself into a FLRW geometry. From a 
physical point of view, this scenario raises (at least) two questions. What is the nature 
of the effective cosmological constant during the de Sitter phase of the history of the 
Universe? What is the mechanism which leads to the end of this phase? Among the 
different inflationary models, the first one to be proposed looks specially interesting. 
It consists of a primordial inflationary model describing the birth and grow of the Universe 
as the result of a cooperative mechanism of zero total energy, generated by a quantum 
fluctuation occurring in a flat background. The idea that a quantum fluctuation may be 
the initial seed for the birth of the Universe has first been proposed by Tryon in 1973 
0. It is because the total hamiltonian of gravity coupled to matter is zero (an important 
consequence of the invariance of general relativity with respect to time reparametrisation) 
and the energy contribution of the gravitational conformal modes is negative that it is 
possible to obtain non-trivial semi-classical solutions of Einstein equations where gravity 
and matter emerges from an empty flat space in a cooperative process. 

In the model proposed by Brout et al. matter is described phenomenologically 
through the action of a scalar field non-minimally coupled to gravity. A consistent asymp- 
totic solution of the field equations 0] describes a de Sitter geometry where the role of 
the cosmological constant is played by the (constant) density of the created matter which 
has, due to the energy conservation law, the characteristic equation of state of a cosmo- 
logical constant Pa + Pa = 0. On the other hand, the consistency of the model implies 
that the (effective) cosmological constant A = SvrpA is a free parameter of the obtained 
solutions. It must not necessarily be of order one, what justify a posteriori the semi- 
classical treatment employed. However, the mass of the scalar quanta fields is very large 
1^: Af ~ Q^/^T, this value being almost insensible to the value of the cosmological constant 
P] even when the vacuum polarization effects are taken into account 0; these effects play 
no significative role in the domain of validity of the model (A << 1). On the other hand, 
the high value of the mass of the quanta is the signal that gravitational strong coupling 
must play an important role in the physics of this problem[|lOl. The scalar field used 



in the model must be interpreted as a phenomenological description of complex objects, 
like black holes |jlll . This interpretation, by the way, is confirmed by the analysis of the 
renormalization of the gravitational constant |jlO[ . Indeed, the renormalized gravitational 



constant becomes singular [|T2| exactly for the critical value of the mass field M = Q^/^T, 
what is the origin of the analytical solution of the semi-classical equations which 
describes, for this special value, the transition from a flat space to a singularity. Since the 
cosmological "constant" in this model is described a gas of black holes, their subsequent 
evaporation through the Hawking radiation mechanism may imply a time- decreasing of 
this effective cosmological " constant" , generating hot radiation in the Universe. 

For the moment, there is no consensus about how the inflationary phase ends. How- 
ever, it is possible to distinguish, from the phenomenological point of view, essentially 
two types of transition to an adiabatic era: a first order phase transition, employed in 
the so-called old and new inflationary models, which deals with an effective cosmological 
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constant given by the energy of a scalar field; a second order phase transition, as it be 
discussed here, based on a variable cosmological " constant" . 

The evaporation of the black holes or the infra-red fluctuations of the gravitational 
conformal modes 0] being responsible for the instability of the de Sitter geometry, it 
seems important to extract observational consequences of the model described above 
through the analysis of a simple scenario. In ||15[ analytical solutions of the FLRW type 
for the phenomenological model described above were obtained for any value of the spatial 
curvature {k = 0, ±1). They depend on four physical parameters: the initial value of the 
cosmological constant, Aq, its residual value Aqo, the initial density of radiation in the 
Universe (Tq and the characteristic time of transmutation of the cosmological constant 
into radiation energy r. The choice of a decaying of the cosmological constant into pure 
radiation is motivated in two ways. From the physical point of view, the temperatures 
encountered in this model permit us to ignore the masses of the usual degrees of freedom. 
From the mathematical point of view, the background model is exactly integrable and 
their perturbations easily computable. 

Our purpose in the present work is to study the evolution of scalar perturbations in the 
inflationary model described before. We will restrict ourselves to the case A^o = ctq = 
and k = 0. These hypothesis lead to the advantage that approximate solutions to the per- 
turbed equations can be easily obtained without to oversimplify the physics. Our model 
differs essentially in two ways from the usual inflationary models. First (after having fixed 
the residual gauge freedom allowed by the coordinate transformations), the perturbation 
equations in the framework of this two fluid model (the effective cosmological constant 
and radiation) will result in a third order perturbed equation (analytically solvable in the 
radiation regime and that otherwise can very well be approximated by a 2 + 1 system of 
decoupled equations, also analytically solvable). The second characteristic of the model 
we discuss here concerns the nature of the initial spectrum of the fluctuations. As it will 
be verified later, the fluctuations evolve freely, i.e., they are determined by the Einstein's 
equations, only from the moment where the temperature of the created radiation is signi- 
ficative. Consequently, we will adopt as initial fluctuation spectrum a thermal spectrum, 
instead of a (hypothetical) spectrum inspired by quantum fluctuations of a vacuum state 
(a choice imposed more by metaphysical^ than by physical considerations; see however 



the discussions in |T^, , for instance) 



This paper is organized as follows. In section 2, the field equations and the solution for 
the background are reminded. In section 3 the perturbed quantities are settled out, and 
approximate solutions presented. In section 4, the spectrum of perturbations is discussed. 
Section 5 is devoted to the conclusions. 



■^During the fifties, Bertrand Russell observed: "The accusation of metaphysics has become in philos- 
ophy something like being a security risk in the public service. . . . The only definition I have found that 
fits all cases is: 'a philosophical opinion not held by the present author'." |T9|. 
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2 The field equations 



In order to be complete, let us review briefly the solutions that will be used in our 

perturbative study. The geometry of the Universe is supposed to be of the type flat 
FLRW: 

ds^^-dt'' + a''{t){dx^ + dy^ + dz^) . (1) 

This metric is a reasonable simpliflcation for the primordial stages of the Universe. More- 
over, in our case, due to the existence of an inflationary phase, the curvature can also 
be neglected in later stages of the evolution of the Universe. The unperturbed material 
source of the Einstein's equations is the sum of two comoving fluids: 

Ti: = (? + P) tJ* k+ P K ^ T^fa..). + T(^s.)u (2) 


with P— p{rad.) + Pa, P— Pirad.) +Pa- The terms representing the radiation perfect fluid 
{P{rad.), Pirad.)) are supposed to represent the high energy matter degrees of freedom under 
the form of radiation: 

Pirad.) = ^3^^^ (3) 

where T is the temperature and u is the effective number of degrees of freedom that we 
will take as being those of the standard model {u = 106.75), this value of u constituting, 
of course, a lower bound at the energy considered. The term pA = A(t)/87r = — 
is purely phenomenological. It represents the density of black holes created during the 
cosmogenesis phase and that evaporate later. We will assume it to be of the form : 

A{t) = Aoe-*/^ , (4) 

where the parameter r is the characteristic time of the black hole evaporation which is of 
the order of the hfetime of the created black holes : 

r~^A4«10^ . (5) 

The equation of state of radiation, p(^rad.) — ^P{rad.), permits to obtain a linear second 
order differential equation for the variable Z — {a dot representing a derivative with 
respect to the cosmological time t): 

Z - ^Aoe-*/^Z = . (6) 

The solution of this equation is given by a combination of modifled Bessel functions: 

Z = aKoiwe-^) + /?/o(ti7e-«) , (7) 



where C = ^ and zu — ^yr^Ao. The coefficients a and /3 are given by 

a ^ w{i;,{w) + Io{w)) , (3 ^ -w{K'o{w) + Ko{w)) , (8) 
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and are fixed by the conditions Z{0) = 1 and Z{0) = +y4Ao/3. 

Notice that there was an error of sign in the ref. [^. The scale factor written there 
corresponds to Z < 0, and the plateau appearing in the graph of the entropy as function 
of time is due to the crossing of the singularity Z = 0. The graphs for the entropy and 
temperature as functions of time, for initial conditions corresponding to an expanding 
Universe, are shown on figures 1 and 2. 
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Fig. 1 Logarithmic plot of the ra- 
diation entropy as function of time for 
Ao = 5.510"^ i.e. = 80. 



Fig. 2 Logarithmic plots of the radi- 
ation temperature as function of time 
for An = 5.510"^ i.e. N, = 80. 



3 The perturbed equations 



The perturbed equations for the scalar modes around the background solutions exhibited 
above are obtained through the standard methods. We suppose the metric perturbations 



(9) 



to be synchronous, i.e., such that 9no= 0. The perturbations in density and velocity read 



1 1 

p(t, x) =P {t)+ p {t, x) , u^{t, x) =u^ +u^{t, x) , 



(10) 



with u = (1,0,0,0). All perturbed quantities can be expressed through a Fourier de- 
composition 



9ij (t, x) 
P {t,x) 
u^{t,x) 



1 



^3/2 
1 

1 



TT' 
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hij{t,n)e-'''-''d 



(11) 
(12) 

(13) 



with n.x = Y^in^x\ n being the comoving wave vector of the perturbations. 

Introducing those perturbed quantities into Einstein's equations, we obtain three cou- 
pled perturbed equations: 

[a^hy = a^6p , (14) 
(a^Sp)- + {p + p)a^9 = Snip + p)a% , (15) 

{{p + v)a''ey = -aHp , (16) 
where 9 = n.u and h = —IGira^^J^ihu- The system of equations (|14|,|15|,|16D is underde- 



termined since p = p(rad.) + Pa and consequently 5p = Sp(^rad.) + SpA. We need a specific 
equation for the effective cosmological constant. Hence, we assume 

V;.T(^eL.))«- = -^uiPAun " PaV^u'' = t'^ Pa (17) 

which reduces to u^d^pA = t~^pa, since pa = —pA = A(t)/87r. For adiabatic perturbations 
(<5pA + SpA = 0), it results (r being supposed to be a constant given by the microphysics) 

u'^d^ Pa +dt Pa= r'^ Pa ■ (18) 

Since u' = 0, we obtain Pa= L{x)e~^l'^ . But this solution is a coordinate artifact. Indeed, 
the coordinate conditions g^Q = yet remain satisfied when the residual infinitesimal 
transformations generated by the vector fields 

eo = 0(f) , ei = tpi{x)a^ - di(j){x)a'^ J ^ (19) 

are performed. The choice 0(x) = ^-^^ allows to make 5A(t) = 0. Notice that, in this 



way, the equations (|T4|,|l5|,|T6[) are now gauge invariants with respect to the remaining 
coordinate transformations. Indeed, residual transformations eo = 0, = ipia"^ only 
modify hij{t, x) up to a term whose time dependence is given by a^; hence, h is defined up 
to a constant, leaving h invariant. In other words, in the synchronous gauge go^ = 0, 6pA 



may be fixed to zero and all solutions of the perturbed equations (|T^JT^,|T6|) now acquire 
a physical meaning, in contrast to what happens in the framework of a one-fluid model. 
In this later case, the third order perturbed equation may be reduced, using the residual 
coordinate freedom, to a second order equation. Here, we face a third order equation 
whose three integration constants have physical meaning. 
From now on, we introduce the new set of variables 

T = ip(Tad.)+P(rad.))a^d , TZ = 5 p{^r-ad.) , H = a^fi , (20) 



which permit us to write the system of equations (|T^,|T5|,|T6D as 



^ = ^ , (21) 

Tl = \-87l{p(^rad.) + Pirad.))a^T-(. , (22) 

t ^ . (23) 

6 a 
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The solutions of this system and their cosmological imphcation will be the subject of the 
next section. 



4 Evolution of perturbations 



We are not interested here in the general solution of the system of equation (pl|-p3|) 
governing the evolution of the perturbations, for arbitrary values of their parameters. We 
shall only focus on the values of the parameters that are relevant for the cosmological 
model we consider here. In order that an inflationary model is acceptable, the scale factor 
must be amplified at least by 70 e-folds, before the Universe reaches the GU temperature 
Tgu = 10"'' where radiative era takes place. At this moment, the expansion of the 
Universe is described by the asymptotic expression 

Zas = a-^^a{^-^s) , + 1 , (24) 

Zt It l 

7 ~ 0.577 being the Euler-Mascheroni constant. Using the Einstein's equation Gqo = 
SttToo, we obtain 

iGV-is=[-^^vT'^^] (25) 



1/2 



Hence, a r^\^^^^vT^^^ e^^- ^ 10^2 ^ ^ iq~65 a e-fold number = 75. 

Accordingly, unless for ^ very close to zero, we may neglect the second term in (0). From 
the asymptotical expression of the Kq Bessel function appearing in eq. (||) we may deduce 
the value of vu (i.e. Aq) insuring a given number A'^e of e-fold : 



w = 2Ne + - In ^ + - In 

2 3 2 



1 ol2 2 /Ti4 

2 3 



(26) 



This implies Aq << 1 for a large range of reasonable values of N^, in particular Aq ~ 
1.4 X 10~^ for A'e = 75, what is compatible with the domain of validity of the model, while 
for Aq = l/(1927r^), which corresponds to the assumption of thermal equilibrium between 
de Sitter temperature and Hawking temperature, we obtain N^, ~ 2400. 

According to the range of the values settled by the initial conditions, we may build 
approximate solutions of the perturbed equations. For instance, let us suppose that 71 
scales with a as . From the perturbed equations, using the conservation of the entropy 
of the radiation as the cosmological constant goes to zero, we obtain 7i ~ a^~^ and 
T ~ a^~'/^. From the perturbed equations, it comes out that in this case 7?. ~ ^ ~ a^~^, 
which is negligible compared to ^ p{rad.)a'^'H ~ . Hence, the perturbed equations can 
be approximated by : 

^ = ^ , (27) 

n = ^Pirad.)a^n , (28) 

r ^ ^5 , (29) 

3 a 
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which can be integrated by quadratures. The general solution is 



■H = n.^ + ^j' Z{t')dt' , (30) 



Z{t) Z{t) Ju 

n{t) = z{t)n{t) = ^z{t) - z{t)H{t) , (31) 

T.T,,'^rm,e ,32) 

3 Jti a[t') 



^ = ^^An^ + z{u)ni] (33) 



with 

" " Z{t 

and TZi, Hi, % being the initial values of TZ, TC and T. 

On the other hand, if pa^Ti can be neglected in front of T, the system reduces to 



dn ^ dT -,11 



where C = /j* -^nT\- Its general solution is then 



T = ,cos(^c)+^K.sm(i^() . (36) 

The initial time tj, where the initial values of TZ, T and Ti have to be considered, 
depend on the size of the perturbations. Today the temperature is of the order of To ~ 2 x 
10~^^, the cosmological observational data which are somehow free of astrophysical noises 
mainly concerns sizes Aq from lOOMpc up to 3.000Mpc, i.e., 2 x up to 6 x 

At these scales the Universe can be considered as homogeneous. At the grand unification 
epoch, the size of these fluctuations are shortened by a factor Tq/Tqu = 2 x 10~^^. They 
must still be reduced by a factor a{ti) / aitcu) in order to obtain their values at the initial 
time. Moreover, the evolution of the fluctuations is given by equations (|14|,|15|,|16D only 
once they grow faster than the microscopic interactions can propagate. It is only then that 
their behavior is no more dictated by the microphysics but only by the overall cosmological 
expansion. This implies that the evolution of perturbations with comoving wave numbers 
n and present day size Aq will only be determined by the perturbation equations from an 
initial time Uin) given by d(tj(n)) = n, i.e., 

a{U{,n)) = -T-^e . (38) 

This time ti{n) may be easily estimated by noting that we~^°^ , and thus tue~^" 
(where ^i(n) = tj(n)/2r), depends very slowly on the number of e-folds. Moreover, near 
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ti{n) = equation ( pS]) easily may be linearized. Thus for each value of A^,, it is possible 
to determine a minimal number A^'"*"(A^) of e-folds such that fluctuations of size as 
today emerges from the horizon at t = 0: 



iV"^"(A^) = In 



A* To A^Tq 
m ■ 



+ 1.93 , (39) 



AtttTgu Att'^tTguV^- 
and giving this number, we obtain from eq. (p6D, for a given number of e-folds, 



(40) 



as the time where a fluctuation of size A* has emerged from the horizon. We restrict our- 
selves to positive values of tj(we do not know the physics before t < 0). This implies that 
the minimum number of e-fold between t = and the beginning of the grand unification 
era we will consider is at least 75. 

The decay of the cosmological constant increases very rapidly the temperature (see 
fig. 1), but nevertheless we shall suppose that the density fluctuations are given by the 
laws of the equilibrium statistical mechanics [|TU]: 



The extra factor a~^{ti{n)) is the jacobian that we have to introduce when we consider 
the Fourier transforms with respect to the comoving variables n which are related to the 
physical momentum Uphys at time tj by 



n _2n Tgu 
a{ti) A^ To a{ti) 



^phys \ \ rn ■ l'^^) 



For consistency, we may check that the ratio of the proper energy density of such fluctua- 
tion to the temperature e(n)/Tj(n) is approximatively equals to 1.8 when the fluctuation 
emerges from the horizon, the temperature Tj(n) at this moment being of the order or 
10~^. We may also check that these values are very few sensitive to the values of r and 
N^. Moreover, we also assume the initial spectrum of 6{ti{n)) and h{ti{n)) being given 
by classical physics, i.e., they are not affected by the overall expansion of the Universe. 



Using standard arguments [19, EO], we obtain 



e{%n) ~ vl';'^T-'/\n)n-"^a{U{n))-'''^ , (43) 
h{%n) ~ vl';fT!'\n)n-'''^a{U{n))-''^ . (44) 

Hence, 7l(ti, n) = Z'^(ti{n))6p{ti, n), T{ti, n) ~ nlZ{ti, n) and Tiiti^ n) ~ 7?.(tj, n)/na(ti, n) 
The second relation means that the longitudinal mode h{ti,p) is determined by the gravi- 
tational fluctuations, with a characteristic time of the order of the size of the fluctuations 
considered (of course, one also could ask about the metaphysical content of these consid- 
erations, that we adopt in the absence of a better scenario). 
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With these assumptions, it is possible to verify that at initial times T/a is about 
2 — 5 X 10^ times greater than ^p(rad.)(^^'^ '■ 

a result again almost insensitive to the values of A^^e and r. As a consequence, during a 
first short period of time we may approximate the system of the perturbation equations 
by the eqs ( p7| , p8| , p9D . Hence, one verify that the functions H, IZ and T given by eqs 
(^,^,0) change little in the beginning of the evolution of perturbations. But the ratio 
between the two terms in the right hand side of (|I3p behaves as : 

—p{rad.)a — ~ Svr^ — 46 

and increases very quickly. Hence, for values of ^ very near ^j, the system will be in 
the domain of validity of the first approximate solution (|30|j3T| ). Taking into account the 
uncertainties on the the parameters of the problem (initial spectrum, precise duration 
of validity of the perturbed equations, etc.), we can adopt them as solutions of the per- 
turbed equations. From these considerations, it comes out that at the grand unification 
temperature, the density contrast 5gu = l^pitcu)/ Paul is given by 



pGuZ^u 3 



and all dependence on the wave number n is contained in the value of the integration 
constant p. This constant is given by the sum of two terms (see eq. (|33|), but the first 
one Z{ti)5pi is significantly more important than the second one. Hence, the spectrum 
index depends only on the initial spectrum of density perturbations and on the dates 
where, for a given size Aq, the fluctuations emerge from the horizon. As the temperature 
T{ti,n) is almost constant during the interval of time where the relevant perturbations 
(lOOMpc < Aq < SOOOMpc) cross the horizon while a{ti{n)) n we immediately obtain 
the spectral index 

= 4 + -^ln[(5p(n)/p)2]^2 . (48) 
a Inn 



5 Conclusions 

In this work, we have studied an inflationary model that joins the radiative regime by 
a second order phase transition. This model is inspired in an inflationary primordial 
scenario where the transition from a de Sitter phase to a radiative phase occurs due to 
the evaporation of primordial black holes. It must be considered as a source of a possible 
alternative to the usual inflationary models, and not as a definitive model. Indeed, the 
underlined physics leading to the background geometry is phenomenological and asks for 
more deep fundamental investigations. It offers, however, the conceptual advantage of 
making no appeal to a self-interacting scalar field. The kinetics of the phase transition 
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asks also for a more rigorous argumentation. However, we believe our results not to be 
dependent on the details of the of the model as the spectral index obtained is very little 
sensible to the values of N^, and r. 

Moreover, we have introduced hypothesis on the nature of the initial spectrum of 
the fluctuations. In particular, we employed a thermal spectrum for the fluctuations as 
they emerge from horizon. Notice that, from the conceptual point of view, the spectrum 
more delicate is that of the h but this quantity does not play a significative role in the 
computation of the spectral index, since the gravitational fluctuations is dominated by 
those, much more important, fluctuations of matter, which is governed by TZ. The nature 
of the spectrum oiSp affects the spectral index essentially by the function a(tj) ~ n which 
appears in the jacobian of density fluctuation defining 5pi and on the behavior of Z{ti, n), 
as appears in the expression for the initial conditions. Of course the obtained value of the 
spectral index is not in agreement with observational data. However here we have only 
discussed scalar fluctuations; it remains to see if tensorial fluctuations cannot reconcile 
the model with the expected Harrison-Zeldovitch spectra. But such an analysis reopen 
the question of the choice of the initial spectrum of the gravitational perturbations, a 
question whose answer, in our opinion is beyond our present knowledge of physics. 
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